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Abstract. Two-dimensional simulations of the coarsening process of the 
isotropic/smectic-A phase transition are presented using a high-order Landau-de 
Gennes type free energy model. Defect annihilation laws for smectic disclinations, 
elementary dislocations, and total dislocation content are determined. The 
computed evolution of the orientational correlation length and disclination density 
is found to be in agreement with previous experimental observations showing 
that disclination interactions dominate the coarsening process. The mechanism 
of smectic disclination movement, limited by the absorption and emission of 
elementary dislocations, is found to be facilitated by curvature walls connecting 
interacting disclinations. At intermediate times in the coarsening process, split- 
core dislocation formation and interactions displaying an effective disclination 
quadrupole configuration are observed. This work provides the framework for 
further understanding of the formation and dynamics of the diverse set of 
curvature defects observed in smectic liquid crystals and other layered material 
systems. 
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1. Introduction 

The study of liquid crystalline mesophases has resulted in profound contributions, 
both through innovative technological advancements and through the increase of 
fundamental understanding. Examples of technological applications range from 
display technology to the cutting-edge field of responsive materials [T]. Unlike 
orientationally ordered nematic liquid crystals, higher order mesophases have garnered 
less attention. Lamellar and columnar liquid crystals, with their degree of translational 
ordering, open up an exciting avenue of progress. For instance, the study of columnar 
mesophases could shed light on the early formation of life on Earth [2l. Focusing 
on the simplest of these higher ordered mesophases, the smectic-A mesophase is a 
physical system which serves as a template for striped systems in general. Striped 
systems are present in a diverse number of areas ranging from biological membranes 
to flexible polymer crystallization [3]. 

Much excitement regarding the smectic-A mesophase is due to its sharing of 
symmetries with certain types of self-assembled block copolymer films. This has 
recently been shown by Harrison et al through experimental observation [S] . Given 
the possible applications of these block copolymer systems in photolithography [6l [7] , 
understanding of their kinetics and dynamics is of great importance. Experimental 
work that has been focused in this direction is limited due to the time scales 
(nanoseconds) and length scales (nanometers) involved. These limitations provide 
a relevant application for modelling and simulation approaches, which have been 
demonstrated as a powerful tool for the prediction of phase ordering behaviour in 
similar systems [8]. Applications of these films to enhance current photolithographic 
techniques in the semiconductor industry are key to advancement in this field [9] . 

In addition to industrial motivations, study of the isotropic/smectic-A mesophase 
transition is interesting for contributions to the expansion of fundamental phase 
transition knowledge. Research cross-fertilization has played a key role in theoretical 
advances through the studying of small-scale rapid phase ordering transitions. 
For example, within the frame of defect kinetics and dynamics, experiments and 
simulations of the isotropic/nematic mesophase transition have been used as a 
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template for the study of the formation of the early universe [10 . Additionally, 
a general topological theory of defects has been cultivated through the study of 
liquid crystalline mesophases as a basis |TT]. While defects of these mesomorphic 
phases are well understood from a macroscopic view, the microscopic structure and 
their interaction dynamics are not [TTl [12]. In general, experimental methods are 
currently limited to static macroscopic studies of mesophase defects. While in recent 
years experimental methods have begun to bridge this gap [11 [H [6l [12] , theoretical 
approaches are vital to further study in this field. 

This work, using modelling and simulation in two-dimensions, studies the first- 
order isotropic/smectic-A mesophase transition. Model parameters were determined, 
in part, using experimental data from 12CB (dodecyl-cyanobiphenyl) [13]. Simulation 
of a temperature quench below the lower stability limit of the isotropic phase was 
performed. The resulting coarsening kinetics, defect structure, and dynamics are the 
focus of this work. The simulation results presented for the two-dimensional smectic-A 
system show strong agreement with experimental observations of a system with the 
same symmetries, self-assembled cylindrical block copolymer films |4l[2. In addition 
to experimental agreement, simulation results show a mesoscopically complete view 
of the defect interaction process. As explained in full detail below, we find that 
disclination interactions are facilitated by curvature walls connecting them. These 
curvature walls are also found to be key in the formation of split-core dislocations 
where pairs of interacting disclinations of opposite sign form dislocations of high 
Burger's vector. Quadrupolar configurations of disclinations are found to dominate 
defect interactions in the intermediate coarsening regime. The presence of these 
disclination configurations result in interesting rheological properties of many different 
materials exhibiting smectic ordering [14] . 

This article is organized as follows: a brief background on relevant types of liquid 
crystalline order and defects (Sections 12. 1112. 2p is given, the model and simulation 
approach are explained (Section l2.3p . and simulation results are discussed (Section[3]). 
Two different viewpoints are used to understand the simulation results: i) overall 
defect /texture kinetics and ii) specific defect dynamics. Overall defect/texture kinetics 
addresses multi-body interactions and the overall evolution of the texture of the 
liquid crystalline domain. This type of analysis is most comparable to experimental 
observations of similar systems. The final sections of this article address specific defect 
interactions at time and length scales that are not fully accessible experimentally. 

2. Background and theory 

2.1. Liquid crystalline order 

Liquid crystalline phases or mesophases are materials which exhibit partial 
orientational and/or translational order. They are composed of anisotropic molecules 
which can be disc-like or rod-like in shape. Thermotropic liquid crystals are typically 
pure-component compounds that exhibit mesophase ordering most greatly in response 
to temperature changes. Lyotropic liquid crystals are mixtures of mesogens, possibly 
with a solvent, that most greatly exhibit mesophase behaviour in response to 
concentration changes. Effects of pressure and external fields also influence mesophase 
behaviour. This work focuses the study of a rod-like thermotropic liquid crystals which 
exhibit a first-order liquid crystalline phase transition. 

An unordered liquid, where there is neither orientational nor translational order 
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(apart from an average intermolecular separation distance) of the molecules, is referred 
to as isotropic. Of the many different mesophases observed, the two of interest in 
this work are the nematic and smectic-A. The nematic mesophase exhibits partial 
orientational order, where molecules conform to a local preferred orientational axis. 
Smectic liquid crystals share the orientational order of nematics, but additionally have 
partial translational order. The simplest of the smectics is the smectic-A mesophase, 
which exhibits one-dimensional lamellar order in the direction of the preferred 
molecular orientational axis. It can be thought of as layers of two-dimensional fluids 
stacked upon each other. A schematic representation of these three phases are shown 
in Figure [TJ 




Figure 1. schematics of a) the isotropic, b) nematic, and c) smectic-A 
mesophases. 

Theoretical characterization of mesophase order is accomplished using order 
parameters that adequately capture the physics involved. These order parameters 
typically have an amplitude and phase associated with them. In order to characterize 
the partial orientational order of the nematic phase, a second order symmetric traceless 
tensor can be used: 

Q = S{nn - i/) (1) 

where n is a unit vector along the average orientational axis and 5* is a scalar which 
represents the extent to which the molecules conform to the average orientational axis. 
For simplicity, the expression for the tensor order parameter ([TJ is simplified using the 
assumption of uniaxial orientational order; refer to |15| for full treatment of the biaxial 
case. 

The smectic-A mesophase has one-dimensional translational order in addition 
to the orientational order found in nematics. Characterizing this mesophase can be 
accomplished through the use of primary (orientational) and secondary (translational) 
order parameters together [16| . This is accomplished using the tensor order parameter 
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([T]) and the complex order parameter: 

* = ?Ae'"^ (2) 

where (f) is the phase and -ip is the scalar amplitude of the density modulation. 

The focus of this work is on the isotropic/smectic-A mesophase transition. 
Experimental observations of this transition show that it is first-order, where material 
properties and order parameters change discontinuously across the phase transition. 
Examples of two of the more well characterized materials that exhibit this transition 
are lOCB and 12CB (cyanobiphenyl) [TT]. 




(a) (b) (c) 



Figure 2. Schematics of different layer structures in the vicinity of a) an 
elementary dislocation b) a — ^ disclination (note, the disclination core in the 
bottom configuration is in the vicinity of the elementary dislocation) c) a 
+ i disclination; solid dots correspond to elementary dislocations and ± signs 
correspond to disclination cores of strength i 



2.2. Defects in the two-dimensional smectic-A mesophase 

Defects in mesophase order result in diverse texture and material properties. 
Only the types of defects observed through the two-dimensional simulation of the 
isotropic/smectic-A phase transition will be addressed. These defects include edge 
dislocations and wedge disclinations, shown in Figure [2] These defects can exist 
with different layering structures, where one is more energetically favourable than 
the other. Dislocations have an associated Burger's vector, b, which describes the 
topological charge associated with it. In a right-handed coordinate system, both 
dislocation structures in Figure[2^ have b — +1. Thus the sign of the Burger's vector 
corresponds to whether the dislocation begins (+) or ends (— ) a set of layers, and the 
magnitude corresponds to the number of layers. 

Many experimentally observed smectic defects such as twist dislocations and the 
diverse set of focal conic domains cannot be captured using two-dimensional simulation 
of a full order parameter model. One type of composite defect is observed, the split- 
core dislocation or giant dislocation [T8| , shown in Figure |3^. This type of dislocation 
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typically has a Burger's vector greater than one, where the dislocation core is split into 
a pair of ±i disclinations. A single edge dislocation of Burger's vector greater than one 
is energetically unstable given the resistance of the layers to expansion/compression 
|19] . Although, as shown in Figure an array of elementary dislocations can result 
in a topologically equivalent defect structure. 




(a) (b) 

Figure 3. a) illustration of a lone it^ disclination pair interaction that forms a 
split-core dislocation where plus and minus signs correspond to ±i disclinations, 
solid circles correspond to elementary dislocations, and hollow circles correspond 
to dislocation content of an attached curvature wall (based on Figure 2a from 
[20) ) b) illustration of topologically equivalent array of elementary dislocations ; 
b = +5 for both of these defect structures. 

Referring to Figure[3^, a curvature wall is attached to the interacting disclinations 
and relieves dilation/compression forces resulting from the disclination cores [18[ I20j . 
Curvature walls separate two disoriented smectic-A domains where the disorientation 
of the two domains is weak [18]. They belong to a class of symmetric planar wall 
structures observed in the smectic-A mesophase which are outlined in [18] . 

Little progress has been made in understanding the role of curvature walls in 
the dynamics of interacting smectic disclinations at the mesoscopic level. Considering 
their key roles in the composition of the rich set of curvature defects resulting from the 
lamellar order of smectics [18] , there is motivation for studies focusing on the dynamics 
of these processes at the length and time scales at which they occur. Theoretical 
approaches, in particular [20], have provided insight into the general structure of these 
wall defects. Experimental observations have, at a macroscopic level, provided insight 
into the role of these wall defects in the coarsening process of similar lamellar systems 
[211 [22] and more recently their role in FCD defects in smectics ^23^. In this context, 
an approach based upon modelling and simulation is appropriate for accessing these 
phenomena at the length and time scales of interest. 

2.3. Model and simulation approach 

A Landau-de Gennes type model for the first order isotropic/smectic-A phase 
transition is used that was initially presented by Mukherjee, Pleiner, and Brand [24[[25] 
and later extended by adding nematic elastic terms [26l [271 EH] : 

fb^la{Q:Q)-h{Q-Q):Q + {Q : Qf + a + /? |vl/|4 
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- ij(g:Q) 1*1' - ieQ : (VWvf) 

a= ao (T - T/v/) ;a = ao (T - Ta/) (3) 

where terms 1-5 are the bulk contributions to the free energy, terms 6-7 are coupUngs 
of nematic and smectic order; both the bulk order and coupling of the nematic director 
and smectic density-wave vector, respectively. Terms 8-10 are the nematic and smectic 
elastic contributions to the free energy, respectively. The order parameters are defined 
in (IT])-([2]), T is temperature, Tni/Tai are the hypothetical second order transition 
temperatures for nematic/smectic-A phase transitions (see |17|). and the remaining 
constants are phenomenological parameters. 

There are four aspects of this model which are key in the characterization of 
smectic-A order: i) use of full tensorial and complex order parameters, ii) inclusion 
of a second-order derivative of the complex order parameter, iii) inclusion of a 
nematic elastic term, and iv) imposition of the layer spacing through coupling with 
orientational order (term 7). Clearly, the use of the full order parameters ([l])-© 
characterising orientational and translational order capture both continuum and phase 
transition phenomena. Additionally, term 10 in the free energy density ([3]) directly 
penalizes smectic layer curvature and is required both for stability and to capture 
the layer curvature contribution where the director deviates from the density-wave 
vector. Preliminary simulations of two-dimensional droplet growth have shown that 
nematic elastic terms are also required due to nematic pre-ordering of the smectic 
front and to properly resolve disclination and elementary dislocation cores. Past 
work has shown that inclusion of the h and I2 terms are necessary to resolve the 
isotropic/nematic interface, but we neglect the I2 term in that our initial droplet 
configurations preclude the possibility to grow to length scales at which heterogeneity 
in the nematic/isotropic interface occurs [29j . Finally, by using a formulation where the 
layer spacing is implicitly determined through the competition of coupling and elastic 
terms, simulation results show that this allows for a strong correlation to experimental 
observations (see Figure H]). 

Material dependent coefficients are also an integral part of modelling and 
simulation using a phenomenological-type model. The coefficients used in this 
work were determined by incorporating experimental data and the use of nonlinear 
programming, see [13] for more details. The use of physically realistic model 
parameters is clearly beneficial for comparison with experimental observations. 

The Landau-Ginzburg time-dependent formulation is used to capture the kinetics 
of the phase transition. Due to the higher order derivative term in the free energy 
functional a higher order functional derivative must be used. The general form of the 
time-dependent formulation is as follows : 

dv 5fb dfb (djb\ „„ / djb \ 



where /i is viscosity, 77 is the order parameter, and fb is the bulk free energy. 

A square computational domain of 3.8 x 10~^^m^ (approximately 50 layers) 
with full periodic boundary conditions was used to simulate bulk conditions. The 
initial conditions for the simulation were a set of smectic droplets of equal diameter 
(approximately one smectic layer, Anm) randomly positioned and oriented throughout 
the domain. The initial value for S, ip, and the layer spacing were set to the 
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Figure 4. Sample computational domain where the surface corresponds to 
-Re(^), the dimensionless density modulation, the contour lines are guides for 
the eye to show defects in smectic order, and periodic boundary conditions 
are shown by the fine vertical centreline; the selections outline different types 
of observed defects (see Section 12.21 for more details): 1) a — ^ disclination 
2) a +i disclination 3) an elementary dislocation 4) an angular quadrupole 
configuration of disclinations (see Section I3.4|l 5) a split-core dislocation 6) 
array of elementary dislocations; the phenomenological coefficients used in this 
simulation were T = SSO/C, T]v/ = 322.85A', Tai = 330. 5X, ao 



: 2.823x10 



7 J 



b = 2.823x10'^, c= 1.972 xlO^^r, qq - ..^^^^.^ , 
8 = 9.792 X lO'';^, e = 1.938 x lO'^pAf, «i = 1 x lO-^^^^^ fej = 1 x lO'^^pAf, 
62 = 3.334 x 10"^*^ Jm, and the ratio of the rotational and diffusional viscosities 
= 25. The layer spacing is approximately inm which applies to all figures 
presenting simulation results in this work. 



1.903 Xl0"5. 



J 



2 X lO^^fTv, 



■■ 3.956x10*^ 
-12, 



homogeneous values for the corresponding temperature of 330if . The initial droplet 
number density was set at 3.287 x 10^^m~^ and the random positioning algorithm 
was constrained to avoid initial impingement. A sample of the simulation results is 
shown in Figure |31 

3. Results and discussion 

Figure [4] shows a sample texture at intermediate times in the coarsening process 
where the types of defects observed include elementary dislocations, ±i disclinations, 
curvature walls (of weak disorientation), and tilt walls (of strong disorientation) 
[TSl[5ni[3Tl[T^. Elementary dislocations, with a Burger's vector magnitude equalling 1, 
are seen throughout the computational domain. An important topological observation 
is that these smectic defects do not solely account for the total dislocation content 
of the domain. Split-core dislocations, formed from pairs of smectic disclinations of 
opposite sign, account for the majority of the dislocation content. This preference for 
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split-core versus elementary dislocations is due to the energy associated with each, 
scaling proportional to /o5(|6|) and 6^ [T^, respectively, where b is the Burger's vector 
of the effective dislocation. 

Further study of Figure 2] shows that curvature walls play an important role in the 
interactions of disclinations, and subsequently, split-core dislocations. These split-core 
dislocations are composite defects and their dynamics are diverse. The results of this 
work are presented sequentially, beginning with base disclination phenomena, followed 
by defect/texture kinetics in the coarsening process, and finally, with the dynamics of 
split-core dislocations. 

3.1. Disclination formation, movement, and basic pair annihilation 

The initial condition of randomly positioned and oriented droplets results in the 
formation of disclinations via the Kibble mechanism [32] . The impingement of 
disoriented domains can be resolved in two general ways in a two-dimensional smectic- 
A system: disclination/curvature wall formation or tilt wall formation (arrays of 
elementary dislocations). The formation of the disclination/curvature wall structure 
is energetically favourable, which is shown by their frequent occurrence as seen in 
Figure [H 

Following the growth period, defects begin to interact and are set in motion 
as a result of these forces. The motion of elementary dislocations is well studied 
|33[ [34], but the dynamics of smectic disclinations are understood only from a 
topological point of view. Disclination motion is observed to require the absorption 
and emission of elementary dislocations. The role of attached curvature walls, which 
relieve layer dilation/compression in the vicinity of the disclination core, will be 
discussed later in Section [575] Referring to Figure [5l the annihilation of two oppositely 
charged disclinations is limited by the topological constraints imposed by the smectic 
layers. The movement of each disclination involves the absorption or emission of an 
elementary dislocation and a transition of the core of the disclination. This process 
is seen for both types of disclinations in Figure \5\ consistent with the predictions 
of topology theory and confirming the general coarsening mechanism proposed by 
Harrison et al 5 . 

The annihilation of a disclination pair results in the emission of similarly charged 
elementary dislocations (resulting from disclination movement). This relationship 
implies that the annihilation rates of each set of defects should differ, with that of 
elementary dislocations being the slower of the two. Additionally, the total dislocation 
content of the domain must include the sum of elementary dislocations and dislocations 
associated with curvature walls (Figure [3^). This dislocation content, which is directly 
related to the existence and interaction of disclinations, is accounted for by including 
defects in smectic order contained by curvature walls (see Figure [3^) . Following the 
disclination annihilation shown in the third frame of Figure [5l the resulting elementary 
dislocations emitted from the disclination annihilation process are free to interact. 
This process is shown in the last three frames of Figure [H] where one of the emitted 
elementary dislocations annihilates with an oppositely charged neighbour. 

3.2. Defect/texture kinetics 

Following the growth regime, multi-body interactions of defects dominate the texture 
dynamics. The approximate self-similarity of smectic texture evolution in time, an 
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Figure 5. The annihilation of a +^/ — ^ disclination pair where: i) t = 
(relative), the initial configuration where the effective Burger's vector of the local 
area is 6 = —1 and the +i disclination has just emitted a —1 dislocation ii) 
t = 500, the — ^ disclination absorbs a +1 dislocation (above it) resulting a 
decrease of the separation distance of the pair by one layer iii) t = 750, another +1 
dislocation is absorbed by the — i disclination decreasing the separation distance 
of the pair by one layer iv) t = 1000, the disclination pair annihilates emitting a —1 
dislocation v) t = 1500, the emitted dislocation annihilates with a +1 dislocation 
in its glide plane vi) t = 2000, the post-annihilation region conserves the overall 
Burger's vector of b = —1; black/white dots are guides for the eye to show + \/^\ 
disclinations. 

inherent characteristic of chaotic natural systems, can be characterized by a governing 
power law of the form t". Due to the periodic nature of the computational domain, an 
infinite smectic continuum is only approximated. As the characteristic length of the 
texture approaches the periodic length scale, the system is no longer self-organized and 
deviates from scale-independence. In this intermediate regime, the system organizes 
into a meta-stable state where the resulting texture evolution stagnates. 

In order to accurately characterize the coarsening process, the orientational 
correlation length was determined directly from the tensor order parameter ([3]) [35j. 
The growth, coarsening, and meta-stable regimes were determined from this result, 
shown in Figure [5^. In addition, the texture length scale (the average separation 
distance between disclinations) was calculated using the disclination density, p~ 2 , and 
both quantities were found to exhibit power law behaviour with n = j. This implies 
that the coarsening process is dominated by disclination interactions. Furthermore, 
this shows that the disclination interactions deviate from those observed in the nematic 
phase which yield an exponent of n = ^ for the texture length scale evolution 
[35] . These results are in good agreement with previous experimental work [4l [5] 
and the deviation from the nematic phase is attributed the topological constraints of 
disclination movements and interactions, which was shown in Section f3.il 

The density of disclinations, elementary dislocations, and total dislocation 
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dimensionless time dimensionless time 

(a) Cb) 



Figure 6. a) (bottom) Orientational correlation length and (top) texture length 
scale (p""'^, where p is the total density of disclinations) versus time; the fit curves 
are f"-^*^ and t'^'^'^, respectively b) defect densities versus time for disclinations 
(bottom), elementary dislocations (middle), and total dislocation content (top); 
the fit curves are t~^-'^^, t~^-^^. and t~^-^^, respectively. The dotted lines in both 
subfigures denote the transition away from the self-similar coarsening regime to 
what is referred to as the intermediate coarsening regime. 



content was determined and is shown in Figure [5)3. The disclination density is 
found to be governed by t~2 power law, agreeing with experimental observations 
[HIS]. Additionally, both dislocation densities are found to obey a t~3 power law, 
differing from that of the disclinations. This difference is qualitatively intuitive 
from observations in the previous section, that annihilation of disclinations can occur 
without changing the net dislocation content. This result is also supported by past 
experimental observations which found that dislocations annihilate at a slower rate 
than disclinations [5]. 

Past modelling and simulation studies of other lamellar pattern-forming systems 
have shown that dislocation annihilation rates obey a power law where n — 
[55] . The effects of the incorporation of fluctuations were also studied showing that 
they play a role at intermediate time scales in the coarsening process and at quench 
depths approaching the bulk transition temperature [36j . Thus the incorporation 
fluctuations would not have a substantial effect on the coarsening results found here 
in that the quench depth is below the lower stability limit of the isotropic phase and 
the orientational correlation length was used to determine the self-similar coarsening 
regime (Figure [6^). 

In the present model, couplings between the phase ordering dynamics and 
momentum balance equation are also neglected. This excludes convection effects 
due to defect-defect interactions and re-orientation induced flow. For example, it is 
known that interaction between two disclinations creates asymmetric vortex flows that 
renormalize the rotational viscosity [37] . Nevertheless, in the present work which is 
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Figure 7. The annihilation of a pair of oppositely charged disclinations facilitated 
by the presence of a curvature wall; (i) t = (relative), the initial configuration 
where the +i disclination begins the emission process an elementary dislocation 
into an attached curvature wall (ii) t = 5000, the +i disclination approaches 
a nearby — i disclination through the absorption and emission of elementary 
dislocations into an attached curvature wall, forming a split-core dislocation 
structure (iii) t = 10000, the disclination pair annihilates emitting their total 
dislocation content into an attached curvature wall; black/white dots are guides 
for the eye to show +^/— § disclinations and contour lines outline defects in the 
smectic order. 

focused on texturing statistics of smectic-A materials, the degree of agreement between 
the presented simulation results and experimental observations strongly suggests that 
neglecting convection does not result in deviation from actual texturing statistics. 

3. 3. Role of curvature walls in disclination interaction 

At intermediate times (see Figure |6K) the coarsening process deviates from the 
previously determined power law behaviour as the texture length scale approaches 
the imposed periodic length. At this point, the number of interacting defects becomes 
small and basic interactions of smectic disclinations are more easily observed. These 
interactions provide insight into the structure of many physical systems that related 
to smectic-A and other lamellar ordered systems. 

The texture structures that are found to form at intermediate times clearly show 
the importance of curvature walls to the movement disclinations and elementary 
dislocations. Disclinations of opposite signs tend to pair and form split-core 
dislocations of high Burger's vector, shown in Figure [T] Interactions of these split- 
core dislocations, which are essentially quadrupolar disclination interactions, have 
been observed experimentally in other mesomorphic systems including cholesterics 
[T5] . nematics under an applied field [38^, liquid crystalline polymer melts [T3], and 
most recently cylindrical block-copolymer films [H [5]. These recent experimental 
observations illuminate the importance of understanding these complex interactions. 

Attached curvature walls play a key role in the dynamics of split-core dislocations. 
Figure [7] shows the formation of an unstable split-core dislocation which eventually 
disassociates into an array of elementary dislocations. The role of attached curvature 
walls as a low energy pathway for the movement, emission, and absorption of 
elementary dislocations is clearly shown. Referring to Figure [71 the curvature wall on 
the right side of the +^ disclination begins with no net Burger's vector and transitions 
to a state with 6 > 0. The curvature wall to the left acts as a low-energy pathway for 
elementary dislocation movement to the -|-^ disclination. Through the absorption and 
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emission process, the +^ disclination approaches a neighbouring oppositely charged 
dischnation and the pair annihilates. In this case the split-core dislocation structure 
was transient and unstable. 

3.4- Quadrupolar interactions of disclinations 

Interactions between disclinations are complex due to the topological constraints 
imposed on their movements by lamellar order. Experimental observations have shown 
that these interactions are essentially quadrupolar interactions of disclinations [5j. 
Utilising the approach by Hudson and Thomas for the characterization of quadrupolar 
disclination interactions in an applied field [38 ' , different configurations of quadrupolar 
interactions are shown in Figure [8] The angular quadrupole is the lowest energy 
formation in that it results in the greatest amount of shielding of like-signed defects 
[38] . This topological configuration is also favourable in the smectic-A mesophase 
because it allows for a highly localized disoriented domain which is resolved almost 
entirely by curvature. This explains the preference for this configuration observed 
in Figure [H The majority of disclination interactions at the intermediate time-scale 
occur in an angular quadrupole configuration. The idealized version of this structure 
would be composed of split-core dislocations of equal Burger's vector magnitude and 
would annihilate isomorphically. 




(a) (b) (c) 



Figure 8. Schematic of the three different configurations of quadrupolar 
interactions: a) effective dipole, b) linear quadrupole, and c) angular quadrupole; 
adapted from 1381 . 

Interactions between split-core dislocations of different magnitudes and in non- 
ideal configurations are more likely. Referring to Figure [H all of the angular 
quadrupole configurations are composed of non-equal pairs of split-core dislocations. 
When such a mismatch occurs, a fully isomorphic annihilation process is not possible. 
Figure [H] shows a disclination quadrupole in a mixed angular quadrupole/effective 
dipole configuration. This interaction shows a mechanism absent from interactions of 
elementary dislocations, where two split-core dislocations partially annihilate. This 
can be thought of as an addition of two split-core dislocations where the inner 
disclination pair annihilates and the by-product elementary dislocations resulting from 
their movement are absorbed by the split-core dislocation that remains. This process 
is found to strongly involve the curvature walls associated with respective split-core 




Figure 9. A pair of split-core dislocations of similar sign but differing Burger's 
vector magnitude (6 = 2 for lower and b = 4 for upper) interact in a mixed angular 
quadrupole/effective dipole configuration (see Figure |8j i) t = (relative), the 
initial configuration where both split-core dislocations are of the same sign (the 
interaction must result in addition/consolidation) ii) t = 5000, the +^ disclination 
in the upper split-core dislocation emits an elementary dislocation allowing glide 
of the pair towards the — ^ disclination of the lower split-core dislocation iii) 
t = 10000, annihilation of the inner pair of oppositely charged disclinations occurs 
joining the respective curvature walls iv) t = 11667, the resulting dislocation 
structure emitted by the annihilated disclinations is absorbed into the curvature 
wall of the remaining split-core dislocation, conserving the total Burger's vector 
{b = 6) v) t = 13333, an initially emitted elementary dislocation climbs toward 
the +2 disclination vi) t = 25000, the additive interaction of the two split-core 
dislocations is complete resulting a single split-core dislocation with Burger's 
vector equalling the sum of the original two; black/white dots are guides for the 
eye to show +5/—^ disclinations and contour lines outline defects in the smectic 
order. 

dislocations. The resulting texture has an equivalent Burger's vector of the sum of 
the two original split-core dislocations. 

The addition mechanism for split-core dislocations is determined from these 
simulation results, but annihilation was not observed. Candidates for such an 
interaction are found to form low-energy angular quadrupole configurations (Figure U]) 
and are meta-stable at the time-scales accessible. Using insight from the results of 
this work and previous experimental observations [5l [38], the mechanism for this 
annihilation process of a pair of split-core dislocations can be deduced. An iterative 
process where single elementary dislocations are emitted and then absorbed by a pair 
of split-core dislocations in an angular quadrupole configuration would result in an 
energy barrier to the process being on the order of the core energy of an elementary 
dislocation. This proposed mechanism is compatible with the constraints observed 
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Figure 10. A decoupled curvature wall self- annihilates: i) t = (relative), 
in its initial state after detaching from a disclination ii) t = 10000, shrinking 
isomorphically iii) t = 20000, after fully annihilating an array of elementary 
dislocations remains; contour lines outline defects in the smectic order and 
periodic boundary conditions are shown by the fine horizontal line. 

both here through simulation and experimentaUy ^5J. This would satisfy the basis for 
the coarsening power law exponent, that disclination interactions are limited by the 
movement of elementary dislocations on the order of the separation distance (of the 
disclinations). 

3.5. Curvature wall instability 

The role of curvature walls in the context of disclination interaction has been discussed, 
but the dynamics of unattached curvature walls has not. The behaviour of this 
wall defect removed of a disclination population to meta-stabilize it is shown in 
Figure 1101 As a pair of oppositely charged disclinations form a split-core dislocation 
(from Figure [7|), one of the curvature walls becomes unattached and isolates a 
disoriented domain. This domain structure is similar to that seen in the angular 
quadrupole configuration of disclinations where the disoriented domain is now resolved 
topologically without disclinations. This domain is unstable in its initial configuration 
and its total area is minimized allowing a local relaxation of the layer curvature. The 
annihilation process is iterative where pairs, along the axis parallel to the smectic 
layers, of dislocations associated with the curvature wall approach and annihilate 
expelling any elementary dislocation content that was previously absorbed. The 
disoriented domain decays in an isomorphic process, conserving the topological charge. 
A meta-stable array of elementary dislocations isolate the remaining area of the 
disoriented domain which is stabilized by a texture-imposed localized dilation of the 
layers. 

This process was not observed in previous experimental work [5] due to the length 
scales at which such interactions occur and the difficulty of identifying transient 
curvature walls. These results suggests that in addition interactions of split-core 
dislocations, isolated curvature wall loops have a role in the texturing process. 

4. Conclusions 

The growth and coarsening process of the two-dimensional first-order isotropic/smectic- 
A mesophase transition was studied, texture growth laws determined, and underlying 
defect interactions explained. The role of smectic disclinations was shown to be key 
in the coarsening process and their dynamics strongly deviating from those in the 
nematic phase. These interactions are summarized: 
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• Topologically constrained smectic disclination interaction is limited by the 
absorption and emission of elementary dislocations (Figure [5]) and this process 
does not necessarily result in a change in total dislocation content (Figure . 

• The coarsening process is governed by the annihilation of disclinations where 
the orientational correlation length obeys a power law with exponent n = 
^ (Figure [S^); defect annihilation power law behaviour for disclinations and 
dislocations differs, where n = for both (Figure [SJd). 

• Curvature walls have a fundamental role in disclination interactions and following 
the multi-body interaction regime, split-core dislocations form from disclination 
pairs of opposite sign and their corresponding curvature walls (Figure [7]) . 

• Base split-core dislocation interactions behave as quadrupolar configurations of 
disclinations which accurately predict their dynamics and meta-stable angular 
quadrupolc configurations (Figures |4] and [9]) . 

• Curvature walls are meta-stable defect structures and if detached from stabilizing 
disclinations self-annihilate in an isomorphic process (Figure llOp . 

The phenomena observed through this modelling and simulation study, such as angular 
quadrupoles of disclinations, may be responsible for many of the unique material 
properties of materials exhibiting lamellar ordering |14j . 

Past theoretical and experimental studies have provided insight into much of the 
phenomena presented, but numerical simulation of finer-grained models is currently 
the only way to study these phenomena further. The simulation results presented 
show that the use of a high-order phenomenological model and experimentally based 
model parameters results in a substantially more complete reproduction of the physical 
smectic-A system. Though current computational resources restrict this initial work to 
two-dimensions, these results show a strong correlation to experimental observations 
of some of the complex curvature defects found in smectic liquid crystals. These 
promising results show that three-dimensional simulation of this model could be used 
to study the formation and dynamics of focal-conic defects on length and time scales 
currently unaccessible via experimental study. 

From the viewpoint of industrial applications, the two-dimensional smectic-A 
system shows a great deal of promise for use in cylindrical block copolymer systems. 
Cutting-edge experimental work on enhancement of the photolithographic process 
using block copolymer films [H [71 [9] could benefit from the presented modelling 
and simulation approach in addition to the generation of material specific model 
parameters [13], incorporation of edge-effects [39], and bulk field terms. 
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